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ABSTRACT 

Two optically obscured Wolf-Rayct (WR) stars have been recently discovered by 
means of their infrared (IR) circumstellar shells, which show signatures of interaction 
with each other. Following the systematics of the WR star catalogues, these stars 
obtain the names WR120bb and WR120bc. In this paper, we present and analyse 
new near-IR, J, H, and K-b'&nd, spectra using the Potsdam Wolf-Rayet (PoWR) 
model atmosphere code. For that purpose, the atomic data base of the code has been 
extended in order to include all significant lines in the near-IR bands. 

The spectra of both stars are classified as WN9h. As their spectra are very similar 
the parameters that we obtained by the spectral analyses hardly differ. Despite their 
late spectral subtype, we found relatively high stellar temperatures of 63 kK. The 
wind composition is dominated by helium, while hydrogen is depleted to 25 per cent 
by mass. 

Because of their location in the Scutum-Centaurus arm, WR 120bb and WR 120bc 
appear highly reddened, Ak b « 2 mag. We adopt a common distance of 5.8 kpc to 
both stars, which complies with the typical absolute -KT-band magnitude for the WN9h 
subtype of —6.5 mag, is consistent with their observed extinction based on comparison 
with other massive stars in the region, and allows for the possibility that their shells are 
interacting with each other. This leads to luminosities of log(L / Lq) = 5.66 and 5.54 
for WR120bb and WR120bc, with large uncertainties due to the adopted distance. 

The values of the luminosities of WR 120bb and WR 120bc imply that the imme- 
diate precursors of both stars were red supcrgiants (RSG). This implies in turn that 
the circumstellar shells associated with WR120bb and WR120bc were formed by in- 
teraction between the WR wind and the dense material shed during the preceding 
RSG phase. 

Key words: line: identification - circumstellar matter - stars: fundamental param- 
eters - stars: massive - stars: Wolf-Rayet. 
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Wolf-Rayet (WR) stars play a key role in the cosmic circuit 
of matter. As they are massive hot stars with a strong stel- 
lar wind, they continuously enrich their environment with 
metals and ionise the surrounding interstellar matter. Fur- 
thermore, they are considered as progenitors of supernovae 
and 7-ray bursts. Nevertheless, the origin of WR stars is 
not yet safely established and their physics is far from being 
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fully understood. Detection and study of new WR stars is 
therefore warranted. 

Until recently, the main channel for detection of WR 
stars was through narrow-band optical surveys (e.g. Shara 
ct al. 1999, and references therein). The advent of infrared 
(IR) observations brought new possibilities. Searches for new 
WR stars using narrow-band mid-IR surveys (Shara et al. 
2009, 2012) and IR colour-based selection of objects (Had- 
field et al. 2007; Mauerhan et al. 2011) led to the almost 
two-fold increase of the known Galactic population of these 
stars. 

Another important tool for revealing WR and other 
types of evolved massive stars is through detection of their 
circumstellar nebulae. Classically it is known that some WR 
stars exhibit ring nebulae that are emitting in Ha and other 
optical lines (Johnson & Hogg 1965; Chu et al. 1983; Dopita 
ct al. 1994; Marston 1995). Because of the high obscuration 
in the Galactic plane (where most sites of massive star for- 
mation are located) the distant WR nebulae can be detected 
mainly through their mid- and far-IR emission (e.g. van Bu- 
ren & McCray 1988; Marston 1991; Mathis et al. 1992), 
whose origin can be attributed to radiatively heated circum- 
stellar dust (van Buren & McCray 1988; Marston 1991). 

Barniske, Oskinova, & Hamann (2008) found a neb- 
ula around the WN9h star WR102ka in the Galactic Cen- 
tre region. They obtained a high-resolution IR spectrum of 
the nebula and analysed the emission from warm dust and 
molecules. To our knowledge, this was the first unambiguous 
detection of a dusty nebula around a WN star (cf. Flagey 
et al. 2011). 

Searches for circumstellar nebulae around evolved mas- 
sive stars with the old generation of IR telescopes (onboard 
the Infrared Astronomical Satellite, the Infrared Space Ob- 
servatory and the Midcourse Space Experiment satellite), led 
to the detection of a number of objects (e.g. Marston 1991; 
Trams et al. 1998; Voors et al. 2000; Egan et al. 2002; Clark 
et al. 2003), some of which were already known from op- 
tical and/or radio observations. Significant progress in de- 
tection of new circumstellar nebulae was achieved with the 
release of the MIPSGAL and other Spitzer Space Telescope 
Legacy Programmes. Using Spitzer data, several groups in- 
dependently discovered many dozens of circumstellar shells 
(Gvaramadze et al. 2010a; Mizuno et al. 2010; Wachter et al. 
2010), while follow-up spectroscopy of the central stars of 
these shells led to the discovery of dozens of WR, luminous 
blue variable, and other massive stars (Gvaramadze et al. 
2009, 2010a,b,c, 2011; Wachter et al. 2010, 2011; Mauerhan 
et al. 2010; Stringfellow et al. 2011, 2012). 

Two WR stars revealed with Spitzer via detection of 
their circumstellar shells are the subject of this paper. 
The WR nature of these stars (formerly known as 2MASS 
J18420630-0348224 and 2MASS J18420827-351039) were 
identified by Mauerhan et al. (2010) by means of IR spec- 
troscopy. In the systematics of the Galactic WR star cat- 
alogue (van der Hucht 2001, 2006) , we call these stars 
WR120bb and WR120bc in the rest of the paper. Mauerhan 
et al. (2010) also identified two more evolved massive stars 
(of spectral types of WC8 and 07— 8 III— I), which along 
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with WR120bb and WR120bc form a previously unknown 
star cluster, provided they lie at similar distances. 

Analyses of WR stars are preferably based on their op- 
tical and ultraviolet spectra, which are not accessible for 
WR120bb and WR120bc because of high interstellar extinc- 
tion. However, attempts to constrain the basic parameters 
of WR stars on basis of their A"-band spectra alone have 
been made (e.g. Barniske et al. 2008; Liermann et al. 2010). 
In the present paper the addition of the J and _ff-bands are 
included along with the K-band analysis, which turn out to 
be very valuable. 

The rest of the paper is organised as follows. In Sec- 
tion 2 we introduce WR120bb and WR120bc and their cir- 
cumstellar nebulae in more detail. Section 3 describes our 
spectroscopic observations. After a short characterisation of 
the Potsdam Wolf-Rayet (PoWR) model atmospheres, the 
spectra are quantitatively analysed in Section 4. The evolu- 
tionary status of the stars and the origin of their nebulae 
are discussed in Section 5. We summarise our findings in 
Section 6. 



2 THE CIRCUMSTELLAR SHELLS MN85 

AND MN86 AND THEIR CENTRAL STARS 

Fig. 1 (left panel) shows the 24 /mi image of two adjacent cir- 
cular shells and their central stars (indicated by circles). The 
shells were discovered in the course of our search for evolved 
massive stars via detection of their circumstellar nebulae in 
the archival data of the Spitzer Legacy Programmes (for mo- 
tivation and the results of this search see Gvaramadze et al. 
2010a). The image was obtained with the Multiband Imag- 
ing Photometer for Spitzer (MIPS; Rieke et al. 2004) within 
the framework of the 24 and 70 Micron Survey of the Inner 
Galactic Disc with MIPS (MIPSGAL; Carey et al. 2009). We 
have called the shells MN85 and MN86 2 (Gvaramadze et al. 
2010a). One of the shells (MN86) can also be seen at 70 /im 
(see fig. 1 in Mauerhan et al. 2010). Both shells were also 
covered by the Galactic Legacy Infrared Mid-Plane Survey 
Extraordinaire (GLIMPSE; Benjamin et al. 2003) carried 
out with the Infrared Array Camera (IRAC; Fazio et al. 
2004) and the Mid-Infrared All Sky Survey carried out with 
the Wide-field Infrared Survey Explorer ( WISE; Wright et al. 
2010). The first survey provides images at 3.6, 4.5, 5.8 and 
8/xm, and the second one at 3.4, 4.6, 12 and 22 /im. The 
shells are invisible in the IRAC wavebands, but both can be 
seen in the WISE 22 /im image (not shown here). There is 
also a hint of 12 /an emission probably associated with the 
nebulae (second left panel of Fig. 1). In Fig. 1 (second panel 
from the right) we also show the WISE 3.4 fim image of the 
same field, where the central stars of MN85 and MN86 are 
indicated by circles and the WC8 and 07-8 III-I stars (see 
Section 1) are marked by diamonds. 

MN85 and MN86 were also covered by the Multi- Array 
Galactic Plane Imaging Survey (MAGPIS; Helfand et al. 
2006) carried out with the Very Large Array (VLA). Fig. 1 
(right panel) shows that the shell MN86 has a clear radio 
counterpart at 20 cm. 

2 In the SIMBAD data base the shells are named [GKF2010] 
MN85 and [GKF2010] MN86. 
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Figure 1. From left to right: images at 24 (im (MIPS), 12 /im {WISE), 3.4/mi ( WISE), and 20 cm (VLA) of the field containing the two 
interacting nebulae MN85 and MN86. The positions of the central (WN9h) stars of the nebulae are marked by circles, while two other 
massive stars (WC8 and 07— 8 III— I) are indicated by diamonds (see text for details). 



The MIPS 24 image shows that both shells are 
brightest where the two shells apparently overlap, which 
might be caused by the interaction between the radially ex- 
panding shells (cf. Gvaramadze et al. 2010a; Wachter et al. 
2010). Indeed, Mauerhan et al. (2010) showed that the over- 
lapping region of the shells is significantly brighter than the 
sum of the surface brightnesses of the neighbouring parts of 
the shells, which indicates that it is intrinsically bright emis- 
sion. The enhanced brightness of this region can also be seen 
in the 70 ^m image of the shells in fig. 1 in Mauerhan et al. 
(2010). If the shells are indeed interacting with each other, 
then their almost circular shape implies that they came in 
contact only recently. The existing data, however, do not 
allow us to unambiguously assert that the physical contact 
between the shells is real. 

The angular radius of the larger shell (MN85) is ~ 
1.9arcmin, while that of MN86 is ~ 1.5arcmin. In Sec- 
tion 4.3 we argue that the distance to the objects is ~ 6kpc, 
so that the linear radius of MN85 and MN86 is ~ 3 pc. Cir- 
cumstellar shells of this size are typical of late WN-type 
(WNL) stars and luminous blue variables. The possible ori- 
gin of MN85 and MN86 are discussed in Section 5. 

The central location of point sources indicated in Fig. 1 
by circles suggests that they might be associated with the 
nebulae. This possibility is supported by the position of 
these stars in the colour-colour diagrams by Hadfield et al. 
(2007), where they fall in the region populated by WR stars. 
The details of the central stars are summarised in Table 1. 
The coordinates and the J, H, K B magnitudes are taken 
from the 2MASS (Two Micron All Sky Survey) All-Sky Cat- 
alog of Point Sources (Cutri et al. 2003). The I magnitudes 
are from the DENIS (Deep Near Infrared Survey of the 
Southern Sky) data base (The DENIS Consortium 2005), 
the IRAC magnitudes are from the GLIMPSE Source Cata- 
logue (I + II + 3D), and the WISE 12 and 22 ^tm magnitudes 
are from the WISE All-Sky Data Release (Cutri et al. 2012). 



Table 1. Details of the central stars associated with the nebulae 
MN85 and MN86. 





WR120bb (MN85) 


WR120bc (MN86) 


RA(J2000) 


18 h 42 m 06?31 


18 h 42 m 08=27 


Dec(J2000) 


-03°48'22'. / 5 


-03°51'02{'9 


I 


28? 4811 


28? 4452 


b 


0?3371 


0?3094 


I (mag) 


17.61±0.14 


16.96±0.10 


J (mag) 


11.95±0.03 


11.85±0.03 


H (mag) 


10.22±0.02 


10.26±0.03 


K B (mag) 


9.16±0.02 


9.27±0.03 


[3.6] (mag) 


8.28±0.04 


8.43±0.04 


[4.5] (mag) 


7.63±0.04 


7.88±0.05 


[5.8] (mag) 


7.43±0.03 


7.67±0.03 


[8.0] (mag) 


7.04±0.03 


7.20±0.03 


[12.0] (mag) 


6.92±0.04 


6.96±0.04 


[22.0] (mag) 


5.23±0.05 


3.85±0.06 



All these catalogues and data bases are accessible through 
the Gator engine 3 . 



3 SPECTROSCOPIC OBSERVATIONS AND 
DATA REDUCTION 

Observations were carried out with the Apache Point Obser- 
vatory (APO) 3.5m telescope using the medium-resolution 
infrared spectrograph TripleSpec (Wilson et al. 2004) , which 
provides coverage from 0.95 to 2.46 /im in five orders with 
a resolving power of about 3200 using a l'.'l x 43" slit. The 
exposures were obtained in a nodded ABBA sequence offset 
along the slit by 20 arcsec to facilitate sky subtraction and 
to allow for the correction of detector artifacts. The slit was 
aligned along the east-west sky projection. WR 120bb was 
observed on 2010 May 19 UT at an airmass of ~1.55 with 

3 http:/ /irsa. ipac.caltech.edu/applications/Gator/ 
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a seeing of ~1 arcsec, and WR 120bc on 2010 May 30 UT 
at an airmass of 1.7 with ~1. 2 arcsec seeing. Peak counts in 
the strongest emission lines for both stars remained below 
~ 14,000 counts, corresponding to a detector non-linearity of 
< 0.1 per cent. 

Flat fielding was performed using dome flat lamps, 
and wavelength calibration utilized the OH sky lines ob- 
tained during the exposures. Sky subtraction, flat fielding, 
wavelength calibration, and extraction of the spectra were 
achieved using a modified version of the software package 
xspextool (originally developed for use with SpeX on the 
NASA Infrared Telescope Facility (IRTF); Cushing, Vacca, 
& Rayner 2004) redesigned for use with the APO Triple- 
Spec specific instrument characteristics. Final wavelength 
calibration in each order resulted in an RMS error of < 1.2 
A for WR 120bb, and < 1.5 A for WR 120bc. The AO V star 
HD 171149 was observed immediately before or after observ- 
ing each star for use in performing the telluric and flux cor- 
rections. The difference in airmass between the telluric and 
programme stars were < 0.06. Telluric corrections and ap- 
proximate flux calibrations were performed with the IDL 
routine xtellcor that implements the method developed 
and described by Vacca, Cushing, & Rayner (2003). 

The resulting combined spectrum of WR 120bb has 
a typical signal-to-noise ratio (S/N) in the continuum of 
ss 100, 60, and 4—18 in the K, H, and J-bands, respec- 
tively. Likewise, the values for the combined spectrum for 
WR120bc are 250, 200, and 6-75. The stars are heavily 
reddened, resulting in the much lower S/N at shorter wave- 
lengths. 



4 WR120bb AND WR120bc: TWO ALMOST 
IDENTICAL WN9h STARS 

4.1 Spectral type 

In Figs. 2 and 3 we present the J, H, and if-band spec- 
tra of WR 120bb and WR 120bc. The spectra are very sim- 
ilar to each other and show numerous emission lines of Hi, 
He I and Hell. Hydrogen lines are generally blended with 
the corresponding Hen lines. We also detect emissions of 
N ill and N IV. The absence of carbon and oxygen lines indi- 
cate that both stars belong to the nitrogen (WN) sequence 
of the WR stars. To determine the subtypes of WR 120bb 
and WR 120bc more precisely, we use the equivalent width 
(EW) ratios: EW(Heii 1.012)/EW(He 1 1.083) and EW(Heii 
2.189) /EW(He II + Br 7 ) (see Crowther et al. 2006, and Ta- 
bles therein). Following Crowther et al. (2006), WN stars 
belong to the subtype WN8 if the equivalent width ratio 
EW(Hen 1.012)/EW(Hel 1.083) is between 0.07 and 0.2 
and if EW(Hen 2.189) /EW(He II + Br 7 ) is between 0.1 and 
0.4. If the equivalent width ratios are below these values, the 
stars belong to the subtype WN9. If we compare these ratios 
with the values that we found for WR 120bb and WR 120bc 
(see Table 2), we can assign the two stars to the WN9 sub- 
type. Following the three-dimensional classification for WN 
stars by Smith, Shara, & Moffat (1996), we add a suffix 'h' 
to WN9 to indicate the presence of hydrogen emission lines, 
so that WR 120bb and WR 120bc are WN9h stars. This re- 
sult only slightly differs from that of Mauerhan et al. (2010), 
who assigned both stars to the WN8-9h subtype. 



Table 2. Equivalent widths (EWs; in A) of emission lines and 
their ratios used for spectral classification of WR 120bb and 
WR 120bc. 



EW 




WR 120bb 


WR 120bc 


Hen 1.012 




28.3±7.4 


31.3±2.8 


He I 1.083 




555.5±10.4 


478.1±9.3 


Hell + Bi'7 




122.2±1.0 


122. Oil. 1 


He II 2.189 




3.4±0.3 


5.0±0.2 


EWfHeii 1.012)/EW(Hei 


1.083) 


0.051 


0.065 


EW(Heii 2.189) /EW(He n 


+ Br 7 ) 


0.028 


0.041 



4.2 Spectral analysis and stellar parameters 

To analyse the stellar spectra and to derive the fundamental 
parameters of WR 120bb and WR 120bc, the PoWR models 
for expanding stellar atmospheres are used. These models 
account for complex model atoms including iron-line blan- 
keting in non-LTE (for a detailed description see Grafener 
et al. 2002; Hamann & Grafener 2003, 2004). For the analysis 
of the near-IR spectra, the model atom of H I had to be ex- 
tended compared to the model grid for WNL-stars (Hamann 
& Grafener 2004) in order to reproduce the Brackett-lines 
in the if-band. Other extended model atoms are He I, Hell, 
N in and N IV. The new atomic data are taken from the NIST 
Atomic Spectra Database (Kramida et al. 2012) and from 
the Atomic Line List (van Hoof 1999). 

The normalised emission line spectra of WR stars de- 
pend on two principal parameters: the stellar temperature, 
T t , and the so-called transformed radius, R t . The stellar 
temperature T* denotes the effective temperature, which is 
related to the stellar radius R, and the bolometric luminos- 
ity L via the Stefan-Boltzmann law L — ATrR^aTf, where 
cr is the Stefan-Boltzmann constant and R* is by definition 
the radius where the Rosseland optical depth reaches 20. 
Alternatively, one may quote the radius where the Rosse- 
land mean optical depth of the stellar wind is tr = 2/3. 
Since the winds of our two program stars are very dense, 
this "pseudo-photosphere" is located far out in the wind, 
at about 3.5 72* (see Table 3). Related to that radius, the 
effective temperature T2/3 is much lower than T* (see Table 
3). 

The transformed radius Rt is related to the mass-loss 
rate M and defined by 



Rt — R* 



sfDM 



2500 kms- 



10- 4 M Q yr- 



2/3 



(1) 



introducing the clumping contrast D and the terminal ve- 
locity of the wind Voo . 

In order to achieve a good fit of the electron scattering 
wings, we set D — 10. This is consistent with recent dis- 
cussions of clumping (see Hamann, Feldmeier, & Oskinova 
2008). Note that the effect of D is simply a scaling of the 
empirically derived mass-loss rate, which is proportional to 
D- 1 ' 2 . 

For the velocity field of the wind we adopt the usual 
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Figure 2. Observed J, H, and i\-band spectra (blue-solid line) of WR120bb, compared with the best-fitting model (red-dotted line) 
with the parameters as given in Table 3. For normalisation the observed spectra were divided by the reddened model continuum, and 
then slightly aligned "by eye" for an exact match of the continuum slope. 



/3-law, with the terminal velocity Voa being a free parameter. 
The exponent ft is set to unity throughout this work. 

The basic model parameters T* and Rt are derived from 
fitting the lines of the model spectra to the lines in the nor- 
malised observed spectra (see Figs. 2 and 3). The normali- 
sation is achieved by dividing the flux calibrated observed 
spectra by the theoretical continuum. 

For a systematic determination of those models which 
fit the observations best, we calculated a small grid of models 



in the T, — Rt parameter plane. The terminal wind velocity 
is Woo = 750 kms -1 , as inferred from the P-Cygni like line 
profiles of the He I lines at 1.70 /im (He I 4d 3 D - 3p 3 P) and 
2.06 /xm (Hei 2p 1 P - 2s X S and Hei 7d 3 D - 4p 3 P). The 
determination of the terminal wind velocity is considerably 
facilitated by fitting both H- and i^-band spectra. Now we 
can use two lines with P-Cygni like profiles. As each of both 
lines gives a slightly different range for Voo, the real value 
can be constrained to the overlap of both intervals. 
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Figure 3. Observed J, H, and fC-band spectra (blue-solid line) of WR120bc, compared with the best-fitting model (red-dotted line) as 
in Fig. 2 



As the spectra of WR 120bb and WR 120bc look very 
similar the best-fitting models have nearly identical param- 
eters (Table 3). 

The stellar temperature is adjusted such that the bal- 
ance between the lines from He I versus He n is reproduced. 
The Hell 7-5 line at 1.16 /im in the J-band and the Hell 
9-6 line at 1.48 pm in the //-band are the most important 
indicators for the determination of the temperature. 

This is important to keep in mind, because there are 



other parameter combinations with much lower tempera- 
tures which give well fitting model spectra as well, but only 
models with high temperatures of about 63 kK can repro- 
duce the strengths of these He II lines. 

At the same time, R t is adjusted, which influences the 
strength of the emission lines in general. 

The hydrogen to helium ratio is derived from fitting 
the line strengths of the Brackett series in the //-band. As 
there are few strong hydrogen lines in the /('-band, the anal- 
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Table 3. Stellar parameters for WR120bb and WR120bc. 



Stellar parameters 


WR 120bb 


WR 120bc 


Spectral type 


WN9h 


WN9h 


T» [kK] 


63 ± 2 


63 ± 2 


logiif LRnl 


0.36 ± 0.01 


0.37 ± 0.01 


logL [Lpnl t 


5.66 


5.54 


v^r, \ km s — M 


750 ± 50 


750 ± 50 




10 


10 


R* \Rp)\ 1 

■** L W J 


5.7 


1.9 


loe:M FM^vr -1 ! t 


—4.4 


— 4.5 


log$i [s^ 1 ] 


48.6 


48.6 


R2/3 [R*] 


3.6 


3.4 


T 2 /z [kK] 


33 ± 2 


34 ± 2 


E B -v [mag] 


5.0 


4.7 


Ak s [mag] 


2.0 


1.8 




3.4 


3.1 


M Ks [mag] 


- 6.7 


-6.4 


cf [kpc] (adopted) 


5.8 


5.8 


He [per cent by mass] 


73 ± 5 


73 ± 5 


H [per cent by mass] 


25 ± 5 


25 ± 5 


N [per cent by mass] 


1.5 


1.5 


C [per cent by mass] 


0.01 


0.01 


Fe [per cent by mass] 


0.14 


0.14 



f The error margins of L, R t and M due to the spectral anal- 
ysis are small, but their values scale with the adopted distance d 
which is less certain (see Sect. 4.3). 



ysis of the //-band, which contains most of the Brackett- 
series, was essential to determine the hydrogen abundance. 
For the other elements, there were only weak and no un- 
blended lines, so that we could not determine their abun- 
dances. That is why for other elements mass fractions which 
are typical for Galactic WN stars - N: 1.5, C: 0.01, Fe: 0.14 
per cent (Hamann & Grafener 2004) - are adopted. 

The parameters from the best-fitting model for 
WR 120bb and WR 120bc are compiled in Table 3. The stel- 
lar temperature of 63 kK is unusually high for subtypes as 
late as WN9. As shown in Fig. 4, the domains of the WN 
subtypes in the T, — _R t -plane have a slope. Usually, the WN 
stars with late subtype and hydrogen are located at large R t , 
i.e. their winds are relatively thin (see Hamann, Grafener, 
& Liermann 2006). Our programme stars, however, show 
very thick winds (small R t ). Their location in the T» — R t - 
diagram is consistent with the WN9 subtype derived above 
(Section 4.1) which is based on line ratios. However, with 
the unusually high density in the wind, higher stellar tem- 
peratures are required to reproduce the observed amount of 
ionised helium. 

In addition, one should remember that T* refers by def- 
inition to the radius R, where the radial Rosseland optical 
depth reaches 20. This point of the atmosphere is of course 
not directly seen from outside. The observable photosphere 
is rather located at R2/3 where the Rosseland optical depth 
has dropped to 2/3. The winds of our program stars are so 
thick that R2/3/R* is as large as 3.6 and 3.4, respectively 
(see Table 3). The effective temperature related to that ra- 
dius is only T2/3 = 33 or 34 kK. Thus the location of R, 
depends on the model assumption that the /3-law for the 
velocity continues in the optically thick part of the expand- 
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TVkK 

30 40 50 60 70 80 100 120 150 200 




4.5 4.6 4.7 4.8 4.9 5.0 5.1 5.2 5.3 
log(7VK) 



Figure 4. The positions of WR120bb and WR120bc (overlap- 
ping diamonds) in the log Rt — log T* diagram. The thick lines 
indicate where the classification criteria for selected subtypes 
(labels) match the line ratios that are predicted from theoreti- 
cal models. The underlying model grid was calculated with the 
PoWR code (see Hamann & Grafener 2004, for details). 



ing atmosphere, albeit this cannot be observed, and has no 
influence on the spectral modelling. 



4.3 Luminosity and distance of WR 120bb and 
WR 120bc 

The absolute dimensions of the observed stars cannot be 
determined by spectral analysis only. They are related to 
the distance of the objects which is a priori unknown for 
WR120bb and WR120bc. Therefore, additional assump- 
tions about the distance or the absolute magnitude have 
to be made in order to determine the luminosity, the stellar 
radius and the mass-loss rate. 

The sightline towards WR120bb and WR120bc (I « 
28°) intersects the Sagittarius Arm (located at ~ 1.5 kpc 
from the Sun) and then (at ~ 3 kpc) enters into the Scutum- 
Centaurus (or Scutum-Crux) Arm (e.g. Churchwell et al. 
2009). At a distance of ~ 6 kpc the sightline crosses a re- 
gion of intense star formation (Garzon et al. 1997), which 
is located in the place where the Scutum-Centaurus Arm 
meets the receding tips of the Galactic traditional and long 
bars (Hammersley et al. 2000; Cabrera-Lavers et al. 2008; 
Churchwell et al. 2009). This region spans the Galactic lat- 
itude range from ~ 24° to 32° and contains numerous very 
massive clusters of red supergiants (Figer et al. 2006; Davies 
et al. 2007; Clark et al. 2009; Alexander et al. 2009; Negueru- 
ela et al. 2010, 2011; Gonzalez-Fernandez & Negueruela 
2012) and the giant H II region W43 with a deeply embedded 
(Ay ~ 30 mag) cluster (named in the SIMBAD data base as 
[BDC99] W43 cluster), which contains at least three evolved 
massive stars (the WN7+a/OB? star WR 121a and two O- 
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Figure 5. MSX 8.3 /im image of the Galactic plane centred at I = 28°, 6 = 0°. The positions of three red supergiant clusters (RSGC1,2,3), 
the central cluster of the giant Hn region W43, and the cluster associated with WR 120bb and WR120bc are indicated by circles. The 
Galactic plane is shown by a dashed line. See text for details. 



type giants or supergiants; Blum, Damineli, & Conti 1999) 4 . 
Fig. 5 shows the Midcourse Space Experiment (MSX) satel- 
lite (Price et al. 2001) image of a region of the Galactic plane 
centred at / = 28° ,b = 0° , with the positions of three red 
supergiant star clusters, RSGC1 (Figer et al. 2006), RSGC2 
(Davies et al. 2007), and RSGC3 (Clark et al. 2009; Alexan- 
der et al. 2009), the central cluster of the giant Hn region 
W43 (Blum et al. 1999), and the cluster associated with 
WR120bb and WR120bc (Mauerhan et al. 2010) indicated 
by circles. 

To constrain the distance to WR 120bb and WR 120bc, 
we may adopt plausible values for their absolute magni- 
tude in the if-band. For eleven WN8-9h stars in the Arches 
cluster, located at the known distance of the Galactic cen- 
ter, Martins et al. (2008) obtained a mean K magnitude 
of —6.5 mag with an individual scatter between —6.8 and 
—5.8 mag J However, these WN stars in the Arches cluster 
might be not representative for the stars discussed here, 
since the former are very luminous, but have relatively low 
mass loss rates which reduces their IR excess. 

The calibration of Mk b for Galactic WN stars by 
Crowther et al. (2006) suggests for WN7-9 stars a much 
smaller value of —5.9 mag. This value is the mean for only 



4 The large ionising flux of the W43 star-forming region (~ 10 s 1 
Lyman continuum photons s~ 1 ; Smith, Bicrmann, & Mezger 
1978) is comparable to that of the very massive (~ 10 4 Mq) star 
cluster NGC3603, which suggests that the central cluster in W43 
contains a large number of yet undetected massive stars. 

5 Note that these magnitudes were derived under the assumption 
of a uniform extinction across the Arches cluster of Ax = 2.8 
mag, while it actually ranges from fs 2 to 4 mag (Espinoza, Sel- 
man, & Melnick 2009). Thus, the actual range of M K for WN8-9h 
stars remains unclear. 



three stars, whose parent clusters and associations are be- 
lieved to be known. It is, however, strongly biased by the 
WN8(h)+cc? star WR66, which is either located at a larger 
distance than implied from its possible membership in the 
anonymous association in Circinus (Lundstrom & Stenholm 
1984) , or is a special case of a low-luminosity WN star origi- 
nating from the binary channel. Without this star, the mean 
Mk b for the remaining two (WN7h and WN9h) stars is 
—6.44 mag, which agrees well with our choice. 

Adopting Mk = —6.5 mag and using the i^-band ex- 
tinction derived from the spectral energy distribution (SED) 
fitting (see below), we get distances of d = 5.4 kpc for 
WR 120bb and d = 6.2 kpc for WR 120bc, which imply that 
the stars are located in the Scutum-Centaurus Arm and 
that the radial separation between them is 0.8 kpc. Such 
a large separation, however, is less likely because both stars 
are probably members of a recently discovered star clus- 
ter (Mauerhan et al. 2010) and therefore are most likely 
spatially associated with each other. As Mk of the WN8- 
9h stars in the Arches cluster span quite a wide range, we 
assumed for both stars a common intermediate distance of 
d = 5.8 kpc. The corresponding stellar radii, mass-loss rates, 
luminosities, and absolute magnitudes are included in Ta- 
bles. For the sake of completeness we also give in Table 3 
the hydrogen ionising luminosity, $i. 

Table 3 shows that the stellar parameters we obtained 
are very similar for both stars. However, the reddening to- 
wards WR120bb and WR120bc is quite different. This dif- 
ference can be caused by the non-uniform extinction across 
the cluster, e.g. due to inhomogeneities in the remainder of 
the parent molecular cloud. Alternatively, the enhanced red- 
dening could be caused by a dense clump of the circumstellar 
shell projected along the line of sight towards WR 120bb (cf. 
Gvaramadze 2001). 
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The observed SED was then fitted for each star (con- 
structed from the flux-calibrated spectra and the photomet- 
ric observations from DENIS, 2MASS, IRAC and WISE) 
with the model SEDs (Fig. 6) by adjusting the parameters 
of the reddening curve. For an appropriate treatment of the 
interstellar extinction, we adopted the reddening law from 
Fitzpatrick (1999), which covers the wavelength range up 
to 50 000 A and allows adjustment of two parameters: the 
colour excess Eb-v and the total-to-selective absorption ra- 
tio Rv = Av / Eb-v , which is 3.1 for the diffuse interstellar 
medium but can vary in general between 2.2 and 5.8 (Fitz- 
patrick 1999). For wavelengths above 50 000 A, the redden- 
ing law from Moneti et al. (2001) is used. For the distance 
of d = 5.8 kpc for both stars, good fits to their SEDs are 
achieved with the parameters as given in Table 3. The ex- 
cess seen at 22 /jm obviously reflects the emission from dust, 
since the PSF of the WISE instrument (FWHM = 12") cov- 
ers a part of the circumstellar nebula. 

The adopted distance to WR 120bb and WR 120bc of 
5.8 kpc implies that the parent cluster of these stars is lo- 
cated near the base of the Scutum-Centaurus Arm, i.e. in 
the region of active star formation caused by the interaction 
between the Galactic bar(s) and the Galactic disk (Garzon 
et al. 1997). This inference is supported by a study of the 
distribution of interstellar extinction around I ~ 28° , which 
shows that Ak s slowly increases from 1.2 to 1.6 between 5- 
6 kpc, and then abruptly increases to 2 — 2.5 mag between 6 
and 7kpc (Negueruela et al. 2011, and references therein). 
Note that A Ks towards WR 120bb and WR 120bc is compa- 
rable to that measured in the direction of the red supergiant 
clusters (f» 1.5 ± 0.3 mag). Since it is believed that all these 
clusters are located at a common distance of ~ 6 kpc (e.g. 
Clark et al. 2009; Negueruela et al. 2011), just in front of 
the large extinction wall, it is likely that the parent cluster 
of WR 120bb and WR 120bc is located at the same distance 
as well. 

The actual distance to WR120bb and WR120bc, how- 
ever, could differ from 5.8 kpc. It would be somewhat smaller 
(for the given Mk) if the circumcluster medium significantly 
contributes to the reddening of these stars. In this connec- 
tion we note that Ak s towards the cluster embedded in W43 
(which is located at d ~ 6 kpc) is ~ 3 mag (Blum et al. 1999) . 
Moreover, the absolute JsT-band magnitudes of WR 120bb 
and WR 120bc can differ from those given in Table 3. Assum- 
ing that they are in the same range as those of the WN8-9h 
stars in the Arches cluster (i.e. —6.8 ... — 5.8 mag; Martins 
et al. 2008), one finds that the distance to WR 120bb and 
WR120bc can range from ~ 4.0 to 7.0 kpc, i.e. both stars 
are still located in the Scutum-Centaurus Arm. Correspond- 
ingly, the luminosities of WR 120bb and WR 120bc and their 
R* and M (given in Table 3) can be scaled to different dis- 
tances as oc d 2 , oc d, and oc d 3 ^ 2 , respectively (Schmutz, 
Hamann, & Wessolowski 1989). 



5 PROGENITORS OF WR 120BB AND 

WR 120BC AND THE ORIGIN OF THEIR 
CIRCUMSTELLAR NEBULAE 

The luminosity of WR 120bb and WR 120bc of log(L/L ) w 
5.5 — 5.7 (derived for our preferred distance of 5.8 kpc) im- 
plies that the initial (zero-age main-sequence) masses of 



these stars were < 40 M (e.g. Ekstrom et al. 2012). A 
massive star of this mass evolves through the sequence: MS 
(main sequence) — > RSG (red supergiant) — s* WN, and ends 
its life in a type II supernova. During the MS phase the star 
creates an extended (tens of pc) bubble filled with hot, ten- 
uous gas (e.g. Weaver et al. 1977). Then, during the RSG 
phase it loses a considerable fraction of its initial mass in 
the form of slow, dense wind, which occupies a compact (a 
few pc) region of enhanced density within the MS bubble. 
This region is surrounded by a dense shell, which is created 
by the high-pressure interior of the MS bubble acting on the 
expanding RSG wind, so that the radius of the shell, rrtsG, 
is determined by the balance between the ram pressure of 
the RSG wind and the thermal pressure in the interior of 
the MS bubble (e.g. D'Ercole 1992): 

MrsG^RSG 7 2/5 3/5,-4/5 /„-, 

^4sg ~ (3850^) 2 / 5 MS Po ' { ' 

where Mrsg an d Vrsg are, respectively, the mass-loss rate 
and the wind velocity during the RSG phase, Lms is the 
mechanical luminosity of the stellar wind during the MS 
phase, po — 1.4mH«o, Wh is the mass of a hydrogen atom, 
and no is the number density of the ambient interstellar 
medium. Since WR 120bb and WR 120bc are members of 
the cluster containing (at least) two other evolved massive 
stars (whose progenitors were nearly as massive as those of 
WR120bb and WR120bc), it is likely that all these stars 
contribute to the thermal pressure of the common bubble 
around the cluster. 

Adopting Mrsg = 3 • 10~ 5 M Q yr -1 and ursg = 
15 kms -1 (the figures typical of RSGs), and assuming 
Lms = 3 • 10 36 ergs -1 and no = 1 cm -2 , one has from equa- 
tion (2) that at the beginning of the WR phase (i.e. at t > 4 
Myr) 

„ / Mrsg \ 1/2 I v RSG \ 1/2 
rRSG ~ 3pC (,lO-5M yr-iJ 1 15 kms-J 

( L MS \ " 1/5 /_no_\ ( t \ 2/5 
\3- lO^ergs- 1 ,/ Ucm^J ^4 Myr J 

which is comparable to the linear radii of the shells MN85 
and MN86 (r s h ~ 3pc; see Section 2). 

The subsequent fast WR wind sweeps up the RSG wind 
and creates a new shell, which expands through the region 
occupied by the RSG wind with a constant velocity (e.g. 
Chevalier & Imamura 1983; D'Ercole 1992) 

/ MwR^vrWRSg\ ^ /„n 

V 3A4rsg J 

where Mwr and %r are, respectively, the mass-loss rate 
and the wind velocity during the WR phase. With Mwr 
and wwr from Table 3, one has from equation (3) u s h ~ 
160 kms -1 , which in turn gives the dynamical age of the 
shells of fdyn ~ r'sh/fsh ~ 2 • 10 4 yr. Note that w s h and td yn 
only weakly depend on the distance (both oc <f 1//2 ). The 
small dynamical age of the shells implies that WR 120bb 
and WR 120bc entered the WR phase only recently, which 
is consistent with the observational fact that WRcircum- 
stellar nebulae are associated exclusively with WNL stars 
(Gvaramadze et al. 2009, 2010a, and references therein), i.e. 
with very young WR stars, whose winds still interact with 
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Figure 6. Observed flux distribution of WR120bb and WR120bc (blue line) in absolute units, including the calibrated spectrum and 
the photometric measurements compiled in Tablet, compared to the emergent flux of the model continuum (red/smooth line), in the 
infrared also shown with spectral lines (red/dotted line). The model flux has been reddened and scaled to the distance according to the 
parameters given in Table 3. 



the dense circumstellar medium. On the other hand, the 
small size of the region occupied by the RSG wind implies 
that the circumstellar shells around WR stars are short-lived 
(several tens of thousands of years; e.g. van Marie, Langer, 
& Garcfa-Segura 2005) objects. 

6 SUMMARY 

We analysed new J, H, and 7f-band spectra of two neigh- 
bouring optically obscured Wolf-Rayet stars, WR 120bb and 
WR 120bc, which were revealed via detection of their (ap- 
parently interacting) circumstellar shells with the Spitzer 
Space Telescope and follow-up spectroscopy of central stars 
of the shells. Our analysis of the spectra was based on the use 
of the Potsdam Wolf-Rayet model atmosphere code, whose 
atomic data base has been extended in order to include all 
significant lines in the near-infrared bands. It is shown that 
the use of the J- and ff-band is of great additional value for 
the analysis compared to an analysis which is based on the 
K-bsmd spectrum alone. Despite the late spectral subtype 
of WR 120bb and WR 120bc, we found relatively high stellar 
temperatures of 63 kK for both stars. The stellar wind com- 
position is dominated by helium with 25 per cent of hydro- 
gen. The stellar spectra are significantly reddened, Ak s ~ 2 



mag, which is consistent with the location of the stars in 
the Scutum-Centaurus arm. Adopting a common distance 
of 5.8 kpc, WR120bb and WR 120bc have luminosities of 
log(L/ Lq) = 5.66 and 5.54 and mass-loss rates of 10~ 4,4 
and 10 -4 ' 5 M0yr _1 , respectively. These values have to be 
considered as highly uncertain due to the uncertainty of the 
adopted distance. The inferred luminosities imply that the 
immediate precursors of WR 120bb and WR 120bc were red 
supergiant stars of initial mass of < 40 Mq, which in turn 
implies that the circumstellar shells around both stars orig- 
inate because of interaction between the Wolf-Rayet wind 
and the dense material shed during the preceding red su- 
pergiant phase. 
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